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Abstract 

The He3 system is studied using the adiabatic hyperspherical representation. We adopt the 
current state-of-the-art helium interaction potential including retardation and the nonadditive 
three-body term, to calculate all low-energy properties of the triatomic 4 He system. The bound 
state energies of the 4 He trimer are computed as well as the 4 He+ 4 He2 elastic scattering cross 
sections, the three-body recombination and collision induced dissociation rates at finite tempera- 
tures. We also treat the system that consists of two 4 He and one 3 He atoms, and compute the 
spectrum of the isotopic trimer 4 He2 3 He, the 3 He+ 4 He2 elastic scattering cross sections, the rates 
for three-body recombination and the collision induced dissociation rate at finite temperatures. 
The effects of retardation and the nonadditive three-body term are investigated. Retardation is 
found to be significant in some cases, while the three-body term plays only a minor role for these 
systems. 
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I. INTRODUCTION 



Three-body atomic systems have attracted considerable interest due to the possibility of 
observing an intriguing property generally referred to as "Efimov physics" [if, [2], [J. The 
most dramatic manifestation of Efimov physics is the possibility of an infinity of three-body 
bound states even when none exist for the separate two-body subsystems. This occurs when 
the two-body scattering length a 12 is large compared to the characteristic range Tq of the 
two-body interaction potential, which is typically on the order of tens of Bohr radii. It 
is predicted that, in the limit of large scattering length (therefore zero two-body binding 
energy), the atom-dimer scattering length also becomes large and the three-body system 
acquires an infinite series of bound states (called Efimov states) whose energies form a 
geometric sequence with zero as an accumulation point. These phenomena persist even when 
r o/°-i2 is nonzero, so that we can expect to observe Efimovphysics in a real physical system. 
The theory of Efimov physics was formulated in 1970 l[ O, [J, but was experimentally 
confirmed only in 2006 [4], by exploiting a Feshbach resonance in an ultracold gas of 133 Cs 
atoms to adjust the scattering length. However, the evidence of Efimov physics was seen in 
the three-body recombination rates, not by direct observation of bound states. It is therefore 
clear that understanding three-body recombination is crucial in establishing the extent to 
which Efimov physics can be realized in the laboratory. 

Helium has long been considered to be one of the most promising candidates for seeing 
Efimov physics since the 4 He dimer has a large scattering length (larger than 100 a.u.). 
The theoretical treatment of triatomic 4 He systems is simple compared to other atomic 
species because there exists only one dimer bound state which has zero orbital angular 
momentum 1 = 0. Several sophisticated helium dimer interaction potentials have been 
developed, of which the most widely used one is probably the LM2M2 potential by Aziz and 
Slaman Q]. This potential has been used to calculate various proper ties of the helium dimer 
and trimer jf], Q, S, [9} as well as their scattering observables 10. 11. Q,Q. The HFD-B3- 



FCI1 potential developed by Aziz et al. 
recombination rates of cold helium atoms 
al. 



141 has also 



is 3, 
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:>een used to calculate the three-body 



171 ] . Jeziorska et al. [18j| and Cencek et 



19] have recently developed not only a helium dimer potential, but also the retardation 



correction to the dimer potential and the nonadditive three-body term. Retarded potentials 

nn 

have been used previously [13|, 120] and were found to have clear impacts on the low-energy 
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behavior of helium systems. Experimentally, the 4 He dimer has been observed by Luo et 



al 



21 



22j and by Schollkopf and Toennies 



23 



|. The latter authors were able to see 
not only the helium dimer but also the trimer and tetramer. The 4 He trimer they observed 
was in its ground state, but no experiment has been able to see the excited state. There is 



actually much discussion currentl y w hether one or 



Efimov states, see Refs 
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his work is along the line of the investigations in Refs. 
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35 



32 [33j,l3 
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36j, which dealt with the spectrum of the helium trimer, the atom-dimer elastic 



scattering in triatomic helium systems, and the three-body recombination of cold helium 
atoms. None of these previous works, however, presented all of these quantities. By doing 
so, we hope to benefit those applying effective field theory to this system. Moreover, we 
adopt the current state-of-the-art triatomic helium interaction potential, including the two- 
body term with retardation by Jeziorska et al. jisl . and the nonadditive three-body 
term of Cencek et al. [19!]. Using this complete potential, we compute the bound state 
energies of the 4 He trimer and the cross sections for elastic scattering between 4 He and 
4 He2- We also calculate the rates for three-body recombination 4 He+ 4 He+ 4 He— s> 4 He2+ 4 He 
and collision induced dissociation 4 He2+ 4 He— ^ 4 He+ 4 He+ 4 He at finite temperatures. In our 
previous paper 1^ . we published the recombination and dissociation rates using the HFD- 
B3-FCI1 potential and including the states with total angular momenta from J = to 3. In 
the present work, we extend the calculations up to J = 7 using the new potential mentioned 
above. In addition, we treat the system consisting of two 4 He and one 3 He atoms, and 
calculate both its spectrum and the cross sections for elastic scattering between 4 He2 and 
3 He. We also obtain the rates for three-body recombination 4 He+ 4 He+ 3 He— > 4 He2+ 3 He and 
collision induced dissociation 4 He2+ 3 He— *- 4 He+ 4 He-|- 3 He at finite temperatures for the first 
time. To summarize, we will present calculations of all low-energy properties of triatomic 
helium for the isotopic combinations that include 4 He2. Moreover, we do so using the best 
available potentials and examine the contribution of the usually-neglected retardation and 
purely three-body terms. 

The key ingredient in our numerical calculations is the adiabatic hyperspherical represen- 



tation [6|, [15|, [38]. Enforcing the boson permutation symmetry is simplified using a modified 



version of the Smith- Whitten coordinate system J, |40j. The i?-matrix method [4l|] is 

used to obtain the scattering S-matrix. Using this ^-matrix, we calculate the atom-dimer 
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elastic scattering cross sections 02, as well as the rates for three-body recombination K% 
and collision induced dissociation D%. Note that the rate equation for the density of helium 
atoms in a thermal gas can be written as 

= ~^^ 3n 4 _ ^ K 3 n l n 3 + 2D 3 n 4d n 4 + 2D' 3 n M n 3 , (1) 

where ra 4 , n id and n 3 are the densities of 4 He, 4 He 2 and 3 He, respectively. K 3 
and K' 3 are the rates for the recombination processes 4 He+ 4 He+ 4 He— > 4 He2+ 4 He and 
4 He+ 4 He+ 3 He— ^ 4 He2+ 3 He, respectively, while D3 and D' 3 are the rates for the dissocia- 
tion processes 4 He2+ 4 He— >- 4 He+ 4 He+ 4 He and 4 He2+ 3 He— s> 4 He+ 4 He+ 3 He, respectively. 

This paper is organized as follows. We explain our method and give all necessary formulas 
for calculating the spectra and the collision properties of triatomic helium systems in Sec. [Til 
The results are presented in Sec. IIHI A summary of this work is given in Sec. IIVI We use 
atomic units throughout except where explicitly stated otherwise. 

II. METHOD 

We solve the Schrodinger equation for three interacting helium atoms using the adiabatic 



hyperspherical representation 



15 



]. In the adiabatic hyperspherical representation, we 
calculate eigenfunctions and eigenvalues of the fixed-hyperradius Hamiltonian in order to 
construct a set of coupled radial equations. The bound state energies can be obtained as the 
discrete eigenvalues of these coupled equations. The scattering S'-matrix can be extracted 
from the same coupled equations using the .R-matrix method. The method employed is 



largely the same as detailed in Ref. [15|. We thus give here a brief outline, emphasizing the 
modifications necessary when the three particles are not identical. 

After separation of the center of mass motion, any three-particle system (in the absence 
of an external field) can be described by six coordinates. Three of these can be chosen as 
the Euler angles a, /3, and 7 that specify the orientation of the body-fixed frame relative 
to the space-fixed frame. The remaining three internal coordinates can be represented by a 
hyperradius R and two hyperangles 6 and (p. To define these internal coordinates, we use a 



slightly modified version of the Smith- Whitten hyperspherical coordinates 15|, |39|, |40|, |42|.|43]. 



We first introduce the three-body reduced mass \i and the scale parameters d^ as follows: 



d 2 = {m k /n){m % + m j} 
in, ■ nij ■ ni),- 



where the indices i, j, k are a cyclic permutation of (1,2,3). The mass-scaled Jacobi 



coordinates 



44 



451 ] are defined by 



Pi = (r 2 - fi)/d 12 , (4) 
m\T\ + m 2 r 2 



p2 = d 



12 



r 3 



(5) 



(7) 



mi + m 2 

Here, r*j is the position of the particle i with mass in the lab-fixed frame. We shall write 
the hyperradius R as follows: 

R 2 = p\ + pl Re[0,oc). (6) 

The hyperangles 6 and <£> are defined by 

(pi),, = j Rcos(tt/4 - 6/2) cos(^/2 + <p 12 /2), 
{px) v = Rsm(7c/A - 6/2) sin(y;/2 + <p 12 /2), 
{pi)z = 0, 

= -^008(^/4 - 0/2) sin(^/2 + ^ 12 /2), 
(p 2 )^ = ii!sin(7r/4 - 0/2) cos(<^/2 + <p 12 /2), 

(P2)z = 0, 

where <pi 2 = 2 arctan(m 2 /p). In the case of three identical particles, the hyperangle ip 
can be reduced from [0, 27r] to [0, 27r/3] in which case the interaction potential becomes 
invariant under reflections about <p = tt/3. Then, in this restricted domain, the solutions 
of the Schrodinger equation are automatically either symmetric (bosonic) or antisymmetric 
(fermionic) with respect to exchange of any two particles (this assumes the spin wave function 
is completely symmetric under exchange of any two particles). In the case that two of the 
three particles (we choose particles 2 and 3) are identical, the interaction potential becomes 
invariant under reflections about if = n, and the solutions of the Schrodinger equation are 
automatically either symmetric (bosonic) or antisymmetric (fermionic) in the range [0, 2tt] 
with respect to exchange of the particles 2 and 3. Note that the definition of the h ype rangle 
ip in Eqs. (JTj) is slightly different from that of Lp prev - in our previous publication [15j. The 
relationship between the present and previous definitions is given by <p = ip pvev - — 4ir/3. To 
demonstrate the effect of particle permutations on the hyperangle <p, Fig. [T] shows contour 
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plots of the potential energy surfaces for both three identical particles and two identical 
particles. 

In our hyperspherical coordinates, the interparticle distances are given by 



with 



r 12 = 2~ 1 / 2 d 12 R[l + sin#cos(^ + ^i 2 )] 1/2 , 

r 23 = 2- 1 / 2 d 23 R[l + sin 8 cos if] 1 / 2 , 

r 31 = 2- 1 / 2 d 31 R[l + sin#cos(y? + </?si)] 1/2 , 

(pi2 = 2 arctan(m 2 //i), y?3i = — 2 arctan(ms//x). 



(8) 



(9) 



We rewrite the Schrodinger equation in terms of a rescaled wavefunction, which is related 
to the usual Schrodinger solution \1/ by ip = i? 5 / 2 \I/. The volume element relevant to integrals 



over 



then becomes 2dR sa\2dd8dipda sin /3d/3d'j. The Schrodinger equation for three 



particles interacting through V(R, 8, ip) now takes the form 



1 rf A 2 15 



2n OR 2 2fiR 2 8/ii? 2 



(10) 



where A 2 is the squared "grand angular momentum operator" and its expression is given in 



Refs. 



15 



43]. 



The interaction potential V(R, 8, ip) used in this work is expressed as a sum of three 
two-body terms and a nonadditive three-body term: 



V(R, 8, ip) = v(r l2 ) + v(r 23 ) + v(r 31 ) + w(r 12 , r 23 , r 31 ). 



(11) 



For the helium dimer potential v(r), we use the representation of Jeziorska et at, given in 
Ref. [l^ . The retardation effect can be incorporated as an additional correction term 37] 
to the dimer potential. Due to this correction, the usual van der Waals term Cgr -6 , valid at 
short range, transforms into Crr^ 7 at long range. The 4 He 2 bound state energies E v i = Eqq 
and scattering lengths a 12 calculated with the retarded and unretarded dimer potentials 



are summarized in Table [I] and reproduce those found in Ref. 18| . For comparison, the 
dimer energy with the older LM2M2 potential is —130 mK jg], giving a scattering length of 
IOOA 12]. Comparison with other available potentials shows genera l ag reement in that all 



give a weakly bound dimer, but the energy can vary by a factor of two [12j, |2Cj • There is every 
reason to believe, however, that the current potential is the best available (see discussion in 
Ref. fl8j). No bound state exists for 4 He 3 He or 3 He 2 . 



For the nonadditive three-body term w(ru, ^23, ^31), Cencek et al. [19|] give two repre- 
sentations. One is based on symmetry-adapted perturbation theory (SAPT); the other, on 
supermolecular coupled-cluster th eory with single, double and noniterative triple excitations 
[CCSD(T)]. As mentioned in Ref. 19j, these two representations are very close to each other. 
In test calculations, we find that the three-body problem shows no noticeable difference be- 
tween the SAPT and CCSD(T) representations, and the numerical results agree with each 
other to four significant digits. Consequently, the results we present here will use only the 
SAPT representation of the three-body term. In the remainder of this paper, we will use 
the term "complete interaction" to indicate that retardation and the three-body term have 
been included. 

The first step that must be carried out is the solution of the fixed- R adiabatic eigenvalue 
equation for a given symmetry J n , with II the total parity, to determine the adiabatic 
eigenfunctions (or channel functions) and eigenvalues (or potential curves). The adiabatic 
eigenfunction expansion gives the wavefunction ip(R, SI) [we will write SI = (8, if, a, /3, 7)] in 
terms of the complete, orthonormal set of angular wavefunctions Q v and radial wavefunctions 
F v 

00 

iP(R,n) = Y,F»(R)®.(R^). (12) 

The channel functions $ y are eigenfunctions of the five-dimensional partial differential equa- 
tion 



A 2 



+ 
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+ V(R,9, ( p) 



$ V (R; ST) — U U (R)$ U (R; ST), 



(13) 



2/ii? 2 8/ii? 2 

whose solutions depend parametrically on R. Insertion of ip from Eq. ( fl2l) into the 
Schrodinger equation from Eq. ( 1101) results in a set of coupled ordinary differential equations 



1 



2^ 



2P VV ,{R)^- + Q VV ,(R) 
aR 



F U {R) 

F V ,(R) — EF V {R). 



(14) 



The coupling elements P UU >(R) and Q UU '(R) involve partial first and second derivatives of 
the channel functions Q v with respect to R, and are defined as follows 

d 



and 



P VV ,{R) = U $ U (R;Q) 



Q VV ,{R) = (( $ U (R; fi) 



OR 

d 2 
dR 2 



$v'{R; Si) 



(15) 



(16) 



The double-bracket matrix element signifies that integrations are carried out only over the 
angular coordinates Q. 

When solving the adiabatic equation (fT3"|) . the degrees of freedom corresponding to the 
Euler angles a, 8 and 7 are treated by expanding the channel functions & U {R] fi) on Wigner 
D functions [15j]. The remaining degrees of freedom 9 and <p are dealt with by expanding 
the channel function onto a direct product of fifth order basis splines [46(. We generate 
the basis splines for 9 from 100 mesh points, while we use 80 mesh points for (p. For small 
hyperradii R, a uniform mesh is employed; for large R, the mesh is designed so that it 
become dense around the two-body coalescence points at (9, ip) = (7r/2,7r — ^12), ( 7r /2,vr) 
and (vr/2, 7r — ^31), where the potential surface V(R, 9, <p) changes abruptly. This leads, for 
example, to a total basis size of 34528 in the case of the 3~ symmetry. We autoparallelize the 
computer program in order to run it on multiple CPU cores. Typically, a calculation of the 
32 lowest eigenvalues of the adiabatic equation ({TBI) takes about 5 minutes of wall clock time 
using 8 1.6-GHz-Itanium2 CPU cores on an SGI Altix 4700 supercomputer. We note that 
our approach, unlike many others employed for this problem, does not use a partial wave 
expansion in the Jacobi-coordinate angular momenta. Our results are, in fact, equivalent to 
including a large number of such partial waves. 



The identical particle symmetry was 



Duilt into the adiabatic equations via the boundary 



conditions in (p. For details, see Ref. [15J. For the case of three identical particles, we have 
$„(R; 9,<p=0, a, 8, 7) = ±$ V (R; 9, ip = —, a, 8, 7), for n = ±1, (17) 

and for the case of two identical particles 

<S> v (R;9,<p = 0,a,8,j) = ±®„(R;9,p = 2ir,a,8,j), for II = ±1. (18) 

These conditions ensure that each solution is either symmetric or antisymmetric with respect 
to exchange of any two identical particles (about half of the solutions are symmetric), thus 
eliminating all states of mixed symmetry. The channel functions with appropriate exchange 
symmetry are then extracted in a postsymmetrization procedure as the solutions that satisfy 

I + P12 



2 

for the case of three identical bosons, and 



$ v (R;n))) = l (19) 



$„CR;n))) = l (20) 



for the case of two identical bosons. 

In practice we solve the adiabatic equation ([TBI for a set of about 200 radial grid points 
Ri up to R ~ 2000 a.u. in order to obtain the potential curves U U (R) and the coupling matrix 
elements P uu i(R) and Q vv i(R). For R > 2000 a.u. they are extrapolated using their known 
asymptotic forms. The details of the numerical calculations are explained in Ref. [l5|. 



III. RESULTS AND DISCUSSION 



As predicted by several authors 
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13 



20 



35, 



36 



47 



481 ]. bound states for 4 He 3 



and 4 He2 3 He can exist only for J n = + , and none for J > 0. Since the 4 He dimer has only 
a single I = bound state, atom-dimer scattering, three-body recombination and collision 
induced dissociation are allowed only for the parity-favored cases, that is, II = (— 1) J . In 
the following, we therefore limit ourselves to J n = + , 1~, 2 + , etc... For each of these cases, 
the lowest adiabatic potential curve v — corresponds asymptotically to two atoms bound 
in a dimer with the third atom far away. Since the dimer has only I = 0, the orbital angular 
momentum of the relative motion between the dimer and the atom should be / 12 3 = J- This 
potential asymptotically behaves as 

U (R) - ^Qoo(R) -» £00 + /l ' 23 ^g 2 + fori? -00, (21) 

where Eqq is the energy of 4 He2- All the higher channels v = 1,2,3, ... for each symmetry 
correspond to three-body continuum states i.e., all three atoms far away from each other as 
R — > 00. Recall that in the adiabatic hyperspherical representation the three-body contin- 
uum is rigorously discretized since the adiabatic Hamiltonian depends only on the bounded 
hyperangles. These three-body continuum channel functions converge asymptotically to the 
hyperspherical harmonics. Therefore, the corresponding potential curves behave as 

U V {R)^ A(A ^2 + ^ fori? -00. (22) 

In principle, A can take on any non-negative integer value, but their possible values are 
restricted by the requirements of permutation symmetry. Figure [2] shows the potential 
curves U V {R) as functions of the hyperradius R for 4 He3 and 4 He2 3 He, J n = + . Note 
that we have also calculated potential curves for 4 He3 system with J n = 1~, 7~ , and for 
4 He 2 3 He with J n = 1", 2+ and 3~. 
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Although we will present cross sections and rates as a function of temperature, they 
are not thermally averaged. Rather, the temperature is obtained from the energy by the 
conversion noted in Table [H If needed, the thermal average can be carried out as described 



in Ref . 



491 ] for K 3 or as described in Ref. [50J] for two-body processes 



A. Bound state energies 

The bound state energies can be obtained as the discrete eigenenergies of the coupled 
equations in Eq. f[T4"l) . We solve these equations by expanding F U (R) onto a set of fifth order 
basis splines in R. Typically, a solution of the coupled equations including 10 adiabatic 
channels takes about 1 minute of CPU time using one 1.6-GHz-Itanium2 processor on an 
SGI Altix 4700 supercomputer. 

The "adiabatic approximation" is the solution of Eq. (TT4"|) with only one channel z/, 
which eliminates the sum over v' on the left-hand side of the equation. This leaves a one- 
dimensional Schrodinger equation with an effective hyperradial potential U U (R) — ^Q UU (R) 
that determines the three-body spectrum in the adiabatic approximation. The lowest v = 
energy level obtained by solving Eq. (fT4l) with all nondiagonal coupling elements neglected 
is a variational upper bound to the true ground state energy E . One can also solve Eq. (|T4|) 
neglecting the diagonal coupling term Q VV (R). This corresponds to the hyperspherical equiv- 
alent of the Born-Oppenheimer approximation. The lowest resulting energy for v = gives 
a rigorous lower bound to E . Finally, solving Eq. (|T4"I) with the sums truncated at i/ max 
gives variational approximations to the exact energies. These energies will thus converge to 
the exact energies from above in the limit — > oo. 

As in previous investigations q q y, y, q, q, y, y, y, 0, q. we find tw ° b ° und 

states for 4 He3 and one for 4 He2 3 He. Table [III summarizes these different levels of approx- 
imation to the bound state energies of 4 He3 and 4 He2 3 He as well as the convergence with 
respect to u max . We also include an estimate of the exact energy obtained by exprapolating 
the z/ max -dependent energies using E Vmsix = E CXACt + a/V^ ax . The relatively slow convergence 
in the number of channels is due, we think, to the large repulsive core of the helium inter- 
action potential. To conclude, we have obtained Eq = —130.86 mK and E\ = —2.5882 mK 
for 4 He3 and E Q = —16.237 mK for 4 He2 3 He with the complete interaction potential. For 
comparisons, Nielsen et al. [8( obtained E = —125.2 mK and E\ = —2.269 mK for 4 He3 
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and E = -13.66 mK for 4 He 2 3 He, using the LM2M2 potential. Using the same LM2M2 
potential, Blume et al. (3) obtained Eq = —125 mK and E\ = —2.28 mK for 4 He3. More- 
over, these energies are also in general agreement with those obtained for other two-body 
potentials. 

From the table, it can be seen that more than 30 channels are needed to obtain con- 
vergence to four digits and that 4 He3 converges more quickly than 4 He2 3 He. It is also 
clear that the adiabatic approximation is better than the Born-Oppenheimer as has been 
noted before jf], Q] • Table [TT] also shows the effects of retardation and the three-body terms. 
Comparing the converged results, we see that the complete potential reduces the binding 
compared to the simple pairwise, unretarded potential. In all cases, retardation is more im- 
portant than the three-body term, giving roughly an order of magnitude larger shift in the 
bound state energies. Given the relative importance of retardation, it is clear why the ener- 
gies are shifted upwards since the retardation corrections shortens the tail of the potentials 
from r~ 6 to r~ 7 . Physically, it is also easy to understand why the three-body term is rela- 
tively unimportant once one notes that its leading-order contribution is the Axilrod- Teller, 
dipole-dipole-dipole dispersion term. Because helium is closed shell and tightly bound, its 
polarizability is small, yielding small dipersion terms in general. We conclude, then, that 
the relative importance of retardation and the three-body term found here should not be 
assumed to be a general result. 



B. Atom-dimer elastic scattering cross sections 

The scattering observables were obtained by solving the coupled equations (fT4"l) using 
a combination of the finite element method (FEM) 50j and the i?-matrix method 41]. 
Typically, about 12 adiabatic channels (thus i/ max = 11) are used, and 10 4 elements, in each 
of which fifth order polynomials are used to expand the radial wavefunction, extend from 
R = 5 to 5 x 10 5 a.u. The scattering ^-matrix is then extracted using the i?-matrix method. 
Each energy took less than 1 minute of wall clock time using one 1.6 GHz Itanium2 processor 
on an SGI Altix 4700 supercomputer. We have checked the stability of the S-matrix with 
respect to the final matching distance, number of FEM sectors, and the number of coupled 
channels, and have found our results accurate to two significant digits. 

The cross section for atom-dimer elastic scattering is expressed in terms of the S'-matrix 
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as 



ff 2 = 2^ (T 2 =2^ — £2 l^o-o - 1| • (23) 

j,n j,n K i,23 

Here, is the partial atom-dimer elastic scattering cross section corresponding to the 
J n symmetry, k% >2 3 = [2yUi >23 (i? — Eqq)} 1 / 2 is the atom-dimer wavenumber and ^1,23 is the 
atom-dimer reduced mass /^i j2 3 = v(i\{m 2 + m^)/{mi + m 2 + m 3 ). The event rate constant 
K 2 is related to the cross section by the formula K 2 = ^23^2/1^1,23- 

The atom-dimer scattering length is another useful quantity that has been studied an 



12 



16 



361 ] . It is defined as 



tan5 . . 

a 1)23 = - lim — A 24 

^12,3^0 fc 1)23 

where 5 is the phase-shift for atom-dimer elastic scattering and is related to the diagonal 
^-matrix element by the formula 

S° t = exp(2z5 ). (25) 

Figure [3]^a) presents the partial cross sections a^ 11 for elastic 4 He+ 4 He 2 scattering as 
functions of the collision energy (E — E 00 ) for the J n = + ,l~,...,7~ symmetries. The 
results obtained with the complete helium interaction potential are shown. For comparisons, 
the figure also shows calculations with various combinations of retarded two-body potential 
and three-body term for the J n = + and 2 + symmetries. In the ultracold limit, here 
E - E 00 < 30/iK, a° 2 + is constant and tends towards the value a 2 + — ► Ana\ 23 , with ai j2 3 = 
230 a.u. =120 A. For comparison, Motovilov et al. jl2| obtained cl\ 2 ^ = 115.5A, usine the 
LM2M2 potential. Curiously, the J > partial cross sections decrease as functions of the 
energy for low energy, behaving as a 2 n oc (E — Eqq)^ 1 . From the Wigner threshold law, 
however, one would expect them to behave as a 2 u oc (E—Eqo) j , since J is the orbital angular 
momentum of the atom relative to the dimer. We have found, though, that collision energies 
of some fiK are still too large to recover this threshold behavior and that the cross sections 
reach the threshold regime only at collision energies less than about 1 nK. The effects of 
retardation and the three-body term are almost negligible for J n = + , but for J n = 2 + , 
retardation increases the partial cross section by about 75% for E — Eqq < 10~ 2 mK, while 
the three-body term has a smaller effect. The trend is reversed at collision energies higher 
than about 0.3 mK. 
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In Fig. [3]^b) we present ovf 11 for 3 He+ 4 He2 elastic scattering for the J n = + , 1~, 2 + and 
3~ symmetries. As in Fig. [3](a), we show the results obtained with the complete helium 
interaction potential and those calculated with all combinations of the retarded two-body 
potential and three-body term. For this system, we found 0x23 = 40 a.u. =21 A. For com- 
parison, Sandhas et al. [131 ] obtained di .23 — 21.0A, using the LM2M2 potential. As found 
in the 4 He+ 4 He2 scattering, the collision energies are found to be still too large to show the 
Wigner threshold behavior a^ 11 oc (E — E 00 ) J . The effects of retardation and the three-body 
term are almost negligible. 

Figure H] shows the phase shifts 5o for s-wave atom-dimer elastic scattering as functions of 
the collision energy. The lower curve shows the results for 4 He+ 4 He2; and the upper one, for 
3 He+ 4 He2. These results were obtained from the complete potential, but we did not observe 
any visible difference from the results obtained neglecting the corrections on this scale. 



C. Three-body recombination rates 

With the S'-matrix in hand, we can also calculate the three-body recombination rates 
The event rate constant for three-body recombination for 4 He+ 4 He+ 4 He— » 4 He2+ 4 He is 

■2 



Ks = V = E^/ n = ^E ^u^ I^J 2 - (26) 

Here, K^ n is the partial recombination rate corresponding to the J n symmetry, and 
k = (2fxE) 1 ^ 2 is the hyperradial wavenumber in the incident channel. Sq^ v represents 
scattering from the initial three-body continuum channels [y = 1, 2, i/ max ) to the final 
atom-dimer channel iy = 0) for the J n symmetry. The factor (3!) derives from the number 
of indistinguishable bosonic particles. In the same spirit, the rate for three-body recombi- 
nation 4 He+ 4 He-|- 3 He— > 4 He2+ 3 He is given by 

K, = V = 2! E E ^iV^W l 2 - (27) 
" J,n u=i ^ K 

We have checked the stability of the results with respect to the final matching distance, 
number of FEM sectors, and the number of coupled channels. Typically, we found the 
results to be accurate to three significant digits at 100 mK. 

Figure [5](a) shows the total rate K 3 for three-body recombination 
4 He+ 4 He+ 4 He^ 4 He 2 + 4 He as well as the partial rates K( u for the J n = 0+ 1~ 2+ 7~ 
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symmetries as functions of the collision energy E. The results obtained with the complete 
helium interaction potential are shown while the total rates calculated with the unretarded 
two-body potential and/or the nonadditive three-body term are also presented in the region 
E < 1 niK. At the lower collision energies, K^ n oc E Xmin , where A m ; n is the minimum 
value of A in Eq. 0221) allowed by permutation symmetry. For J n = + , 1~,2 + and 3~ in 
4 He3, we have A m i n = 0, 3, 2 and 3, respectively; for J > 3, A m i n = J. These threshold 
behaviors are predicted by a generalized Wigner threshold law [51]. For collision energies 
E < 30/iK, the total recombination rate K 3 is constant, and the + partial recombination 
rate K® + dominates. In the ultracold limit, we obtain K3 = 9.93 x 1CT 28 cm 6 /s with 
the complete helium interaction potential, while we have 4.45 x 10 -28 cm 6 /s without 
retardation, 1.01 x 10 -27 cm 6 /s without the three-body term, and 4.54 x 10~ 28 cm 6 /s 
without either correction. Hence, retardation increases the recombination rate by about 
120%, while the three-body term has a much smaller effect, lowering the recombination 
rate by about 2%. The present results differ from our previous calculations [jj]] using the 
HFD-B3-FCI1 potential, which gave K 3 = 7.1 x 10~ 28 cm 6 /s. Shepard Q also obtained 
K 3 = 7.09 x 10~ 28 cm 6 /s using the same potential. 

Figure [5(b) shows the total rate K 3 for three-body recombination 
4 He+ 4 He+ 3 He^ 4 He 2 + 3 He as well as the partial rates K^ n for the J n = + ,l",2+ 
and 3~ symmetries as functions of the collision energy E. At the lower collision energies, 
K( n oc E x ™, with A min = 0, 1, 2 and 3, for J n = + , 1", 2+ and 3", respectively [5l|. In the 
ultracold limit, we obtain K 3 = 9.83 x 10~ 27 cm 6 /s with the complete helium interaction 
potential, while we have 8.97 x 10~ 27 cm 6 /s without retardation, 9.78 x 10~ 27 cm 6 /s without 
the three-body term, and 8.98 x 10 -27 cm 6 /s without either correction. Hence, retardation 
increases the recombination rate by about 10%. We observe that the three-body term has a 
smaller effect by far than retardation, lowering the recombination rate by only about 0.05%. 
Interestingly, recombination is much more efficient here with K 3 an order of magnitude 
larger than the homonuclear case above. 
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D. Collision induced dissociation rates 



We can easily calculate the rates for collision induced dissociation since they require the 
same S-matrix elements as recombination. The collision induced dissociation rate is 

D 3 = £»«f = £ Df = £ £ ^^l^„r, (28) 

^1,23 j, n J,n i/=l ^1,23^1,23 

where D^ u is the partial dissociation rate corresponding to the J n symmetry. 

Figure [6](a) shows D 3 for 4 He2+ 4 He— s> 4 He+ 4 He+ 4 He as well as the partial rates D^ n for 
J n = + , 1~, 2 + , 7~. At the lower collision energies, D^ n oc E Xmin+2 , with the same A m i n 



as Ki u |5jJ. Because it uses exactly the same S- matrix elements as K3, the effects of retar- 
dation and the three-body term are similar. In particular, retardation increases -D3 by about 
134% at low energies, but the three-body term has little effect, lowering the dissociation rate 
only by a few percent. In Fig. 0(b), we show D 3 for 4 He2+ 3 He— > 4 He+ 4 He+ 3 He. 

IV. SUMMARY 

In this work, we have carried out a study of triatomic helium systems using the adia- 
batic hyperspherical approach. Adopting the most realistic, state-of-the-art helium interac- 
tion potential available, we update and extend previous investigations on the bound state 
and scattering properties of the 4 He3 and 4 He2 3 He systems. Based on the most accurate 
helium-helium interactions, these three-body potentials also include two-body retardation 
corrections as well as a non-additive three-body contribution. We have thus been able to 
confirm our statement in Ref. pjj] that the three-body term plays only a minor role. In 
contrast, the effects of retardation are significant, increasing the three-body recombination 
rate by a factor of two to three. 

Systems of ground-state helium atoms are relatively simple because there exists only one 
/ = dimer bound states for 4 He. Further computational improvements must be imple- 
mented before we can extend this approach to the more complicated cases of alkali atoms, 
where the larger number of sharp nonadiabatic avoided crossings pose difficulties. Further, 
in alkali systems, we expect that t 



been recognized previously 52, 



:he three-body term will play a much larger role as has 
53]. This work shows, however, that retardation cannot 
be neglected in calculations of ultracold scattering properties based on ab initio potential 
energies, especially for three-body recombination. 
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FIG. 1: (Color online) Contour plot of the potential energy surfaces as functions of the hyperangles 
9 and ip at R = 15 a.u. for a trimer with three identical particles (a) and one with only two 
identical particles (b). In (a), the translation and reflection symmetries of the potential surface 
can be seen at (p = nir/3, (n = 1 — 5), while in (b) we can identify only the reflection symmetry 
at ip = ir. Exaggerated masses are used in (b) to more easily identify the symmetries. The dotted 
line corresponds to the lowest contour line, the dashed line to the second lowest contour line. 
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FIG. 2: (Color online) Adiabatic hyper spherical potential curves U V (R) for (a) 4 He 3 , J n = 0+ and 
(b) 4 He2 3 He, J n = + . The values of A indicate the asymptotic behavior of the potential curves, 
as given in Eq. ([22]) . (a) The 6 lowest potential curves U V (R) for 4 He3. In addition, the two bound 
state energies are shown as horizontal lines, (b) The 9 lowest potential curves U U (R) for 4 He2 3 He. 
The bound state energy is indicated as a horizontal line. 
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FIG. 3: (Color online) Partial cross sections cr^ for (a) 4 He+ 4 He2 elastic scattering in the J n = 
+ ,l~,...,7~ symmetries and for (b) 3 He+ 4 He2 elastic scattering in the J n = + ,l~,2 + and 3~ 
symmetries, as functions of the collision energy (E — Eqq). The results obtained with the complete 
helium interaction potential are shown for collision energies 1(T 3 <(E- E 00 ) < 2 mK, while those 
calculated with various combinations of the retarded two-body potential and three-body term are 
also presented for J n = + and 2 + . Note that the three-body breakup threshold is located at 
E — Eqq = 1.564 mK (with the retarded two-body potential) or 1.728 mK (with the unretarded 
potential). 
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FIG. 4: (Color online) Phase shifts <5o for s-wave atom-dimer scattering as functions of the collision 
energy. The lower curve shows the results for 4 He+ 4 He2 and the upper one represents the results 
for 3 He+ 4 He 2 . Note that the energies are given relative to the two-body dissociation threshold. 
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FIG. 5: (Color online) Total rate K$ and partial rates for three-body recombination, (a) 

4 He+ 4 He+ 4 He— > 4 He2+ 4 He and (b) 4 He+ 4 He+ 3 He— > 4 He2+ 3 He, as functions of the collision energy 
E. In (a), K% and -f^/ n for the J n = + , 1~, 2 + , 7~ symmetries obtained with the complete 
helium interaction potential are shown while K3 calculated with all combinations of the retarded 
two-body potential and three-body term are also presented in the region E < 1 mK. In (b), K% 
and for J n = + , 1~,2 + and 3~ symmetries obtained with the complete helium interaction 
potential are shown along with K3 calculated with all combinations of the retarded two-body 
potential and three-body term. 
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FIG. 6: (Color online) Total rate -D3 and partial rates -D3 11 for collision induced dissociation, 

(a) 4 He 2 + 4 He^ 4 He+ 4 He+ 4 He and (b) 4 He 2 + 3 He^ 4 He+ 4 He+ 3 He, as functions of the collision 
energy E. In (a), D% and for the J n = + , 1 _ , 2 + , 7~ symmetries obtained with the 
complete helium interaction potential are shown while D3 calculated with all combinations of the 
retarded two-body potential and three-body term are also presented in the region E < 1 mK. In 

(b) , -D3 and for the J n = + , 1~, 2 + and 3~ symmetries obtained with the complete helium 
interaction potential are shown along with D3 calculated with all combinations of the retarded 
two-body potential and three-body term. 
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TABLE I: 4 He 2 bound state energy Eqo (1 a.u. = 3.1577465 x 10 8 mK) and scattering length a\2 
calculated with the retarded and unretarded dimer potentials. 



potential 


£ 00 (a.u.) 


E 00 (mK) 


ai2 (a.u.) 


ai2 (A) 


retarded 


-4.953 x 1(T 9 


-1.564 


173.5 


91.81 


unretarded 


-5.472 x 1(T 9 


-1.728 


165.4 


87.53 
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TABLE II: Ground and excited state energies for He3 and ground state energy for 4 He2 3 He at 

various levels of approximation. The v max - ► bound state energies are obtained by extrapolation. 

All energies are given in units of mK, and are relative to the three-body break-up threshold. 

Eq Ei 
retardation yes yes no no yes yes no no 

3-body term yes no yes no yes no yes no 

4 He 3 
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